Many-body interactions and correlations in coarse-grained descriptions of polymer 

solutions 
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We calculate the two, three, four, and five-body (state independent) effective potentials between 
the centers of mass (CM) of self avoiding walk polymers by Monte-Carlo simulations. For full 
overlap, these coarse-grained n-body interactions oscillate in sign as ( — 1)™, and decrease in absolute 
magnitude with increasing n. We find semi-quantitative agreement with a scaling theory, and use 
this to discuss how the coarse-grained free energy converges when expanded to arbitrary order 
in the many-body potentials. We also derive effective density dependent 2-body potentials which 
exactly reproduce the pair-correlations between the CM of the self avoiding walk polymers. The 
density dependence of these pair potentials can be largely understood from the effects of the density 
independent 3-body potential. Triplet correlations between the CM of the polymers are surprisingly 
well, but not exactly, described by our coarse-grained effective pair potential picture. In fact, 
we demonstrate that a pair-potential cannot simultaneously reproduce the two and three body 
correlations in a system with many-body interactions. However, the deviations that do occur in our 
system are very small, and can be explained by the direct influence of 3-body potentials. 
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I. INTRODUCTION 



An efficient statistical description of condensed matter 
systems and materials almost invariably involves some 
degree of coarse-graining, whereby a large fraction of 
the initial microscopic degrees of freedom are traced out, 
leaving a much reduced space of variables associated with 
the composite entities or pseudo-particles. The latter are 
then coupled via effective interactions which result from 
the partial averaging over the initial microscopic degrees 
of freedom. The reduction of the initial multicomponent 
system to a coarse-grained system with a substantially 
smaller set of composite particles implies that the re- 
sulting effective interactions may involve three-body and 
higher-order contributions, even if the original multicom- 
ponent system involved only pair- wise additive forces, 
like Coulombic interactions. Alternatively, one may wish 
to retain the simplicity of pair-wise additivity of the ef- 
fective interactions, but the price to pay is that such ef- 
fective pair potentials are then state-dependent, e.g. are 
functions of the temperature and/or density. The rea- 
son for this is that the effective interaction energy is a 
free energy associated with the averaged-out degrees of 
freedom, which generally has an entropic component. 

There are many examples of the coarse-graining proce- 
dure which has just been outlined. In molecular systems 
the forces between nuclei result from gradients of the elec- 
tronic ground state energy surface which depends para- 
metrically on the nuclear coordinates, and adjusts adia- 



batically to the slow motion of the latter within the Born- 
Oppenheimer approximation. This scenario is mimicked, 
at least at the level of valence electrons, in "ab initio" 
Molecular Dynamics simulations pioneered by Car and 
Parrinello Q| . However in situations where no strong co- 
valent or hydrogen bonding is present, the more common 
route is to represent the total ground state electronic en- 
ergy surface by a sum of one-body, two-body and higher 
order terms. The one-body contribution is the sum of 
the ground state energies associated with individual, iso- 
lated chemical entities (atoms, ions, or molecules). The 
two-body term is made up of the sum of pair potentials 
acting between molecules and higher order terms corre- 
spond to isolated clusters of three or more molecules. 
The sum of ground state energies of individual molecules 
does not contribute to the forces between them, and can 
hence be ignored in the description of collective equilib- 
rium or transport processes which do not involve chem- 
ical reactions. In the simplest case of rare-gas atoms, 
the pair-wise interactions would include overlap repul- 
sion at short range and dispersion forces at long range, 
while triplet interactions would include, among others, 
the Axilrod- Teller triple dipole dispersion potential M, 
which contributes very significantly to the thermody- 
namic and transport properties of the heavier rare gases 
in their condensed states |3||I| . The effect of the higher or- 
der interactions can be approximately incorporated into 
an effective pair potential, which differs from the bare 
pair potential, valid for an isolated pair of molecules, and 
becomes density-dependent f\M- 
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Similarly in metals an effective interaction between 
ions may be determined by tracing out the conduction 
electrons, using perturbation theory or response theory 
fj,||. Treating the ion-electron coupling to lowest order 
(linear response) leads to a structure-independent vol- 
ume term and to a pair-wise screened effective poten- 
tial between "dressed" ions or pseudo-atoms, which both 
depend on the macroscopic conduction electron density. 
The two-body level is generally sufficient for alkali and 
other simple metals, in part because of a quantum inter- 
ference effect which strongly decreases the magnitude of 
the higher order response terms || . For multivalent and 
transition metals many-body effective interactions can no 
longer be neglecte d H , and a full "ab initio" treatment 
may be necessary pif T 

Coarse-graining becomes crucial in the highly asym- 
metric systems of soft matter, involving macromolecules 
or colloidal particles, as well as molecular scale entities, 
like solvent molecules or ions. The latter are traced out to 
derive effective interactions between the electric double- 
layers associated with charged surfaces (colloids, mem- 
branes etc..) frnj |. The microscopic ions play a role simi- 
lar to valence electrons in metals, but with the important 
difference that quantum degeneracy effects are absent, 
and that finite temperature entropic effects control the 
width of the electric double-layers, with a resulting den- 
sity and temperature dependence of the effective interac- 
tions between the mesoscopic colloidal particles. A good 
example of a such a pair potential is provided by the 
classic Derjaguin-Landau-Verwey-Overbeek (DLVO) ef- 
fective pair potential between spherical charged colloidal 
particles . Three-body interactions can be derived 
in similar fashion |l2| . However, the averaging over the 
microion degrees of freedom also leads to a structure- 
independent, but state-dependent volume term, part of 
which is associated with the self-energy of individual 
double-layers, and which is reminiscent of the volume 
term in metals Jl^Jl^] ; this term has a profound effect on 
the phase diagram of charge-stabilized colloids |]l^ , |l5| . 

Another important class of effective interactions of en- 
tropic origin, which follow from averaging over the con- 
figurations of non-adsorbing polymers and small colloidal 
particles, are the depletion forces which have received 
much renewed interest in recent years [[l6 17 1. Deple- 
tion pair potentials depend strongly on the concentra- 
tion of the depletant; recent attempts have been made 
to compute the three-body interactions from simulations 
or density functional theory Volume terms arising 
from the depletion potential picture may have an impor- 
tant effect on the osmotic equation of state, but they 
are not expected to influence the phase behavior in these 
uncharged systems [ fl6f . 

The present paper is concerned with a coarse-grained 
description of dilute and semi-dilute solutions of poly- 
mers in good solvent. Whereas tracing out the micro- 
scopic ions in a charged colloidal suspension has many 
analogies with the liquid metal problem, the coarse- 
graining of neutral polymers, achieved by integrating out 



the internal monomeric degrees of freedom, resembles 
more closely the case of effective potentials between neu- 
tral atoms and molecules, obtained by tracing out the 
internal electronic degrees of freedom. The basic idea 
explored here, which goes back to Flory and Krigbaum 
fll9f , is to represent a set of polymers, each made up of 
L monomers or segments, as single particles, interacting 
with each other through an effective interaction between 
their centers of mass (CM). The important point, real- 
ized by Grosberg et al. |20| , is that the effective pair 
interaction remains finite, even for infinitely long poly- 
mers. Monte Carlo simulations pl| ] and renormalization 
group calculations |2^] show that for two isolated non- 
intersecting polymer coils, the effective potential between 
their CM's is of order 2ksT in the scaling limit, i.e. for L 
going to infinity, while the range of the interaction is of 
the order of the radius of gyration, _R g , of the polymers. 
Recently we have extended this investigation by simulat- 
ing large systems of self-avoiding walk (SAW) polymers 
at finite concentration |23|]24| ]. The resulting CM pair 
distribution function gyr) was then inverted to yield a 
concentration-dependent effective pair potential v(r;p), 
where p is the number of polymer coils per unit volume. 
Although v(r; p) was not found to change dramatically 
with p, the p-dependence of v is very significant for the 
accurate determination of the osmotic properties of di- 
lute and semi-dilute polymer solutions. 

In this paper we adopt a somewhat different point of 
view, by determining state-independent effective pair, 
triplet, quadruplet and quintuplet interactions; these 
many-body interactions are determined by successively 
considering clusters of 2, 3, 4, and 5 SAW polymer chains, 
determining the corresponding n-body distribution func- 
tions, from which an effective n-body potential is derived. 
Triplet interactions between the cores of star-polymers 
have recently been determined in a similar way [^| . The 
next step taken in this paper is to relate the low-density 
(state-independent) pair and triplet interactions to the 
density-dependent effective interactions determined in 
our earlier work p3|p4] |. A final section will be devoted 
to an analysis of three-body correlations as measured by 
bond-angle distributions and deviations from the Kirk- 
wood superposition approximation |2l| . The convergence 
of the series of n-body interactions is assessed on the basis 
of scaling arguments in the Appendix. 



II. SIMULATION MODELS AND METHODS 

Many properties of polymers in a good solvent are well 
described by models which ignore all microscopic details 
of the intermolecular interactions, except their excluded 
volume. For that reason polymers are often modeled as 
self avoiding walks (SAW) on a lattice [p7|-p9[, a model 
lending itself well to efficient computer simulations. We 
consider the situation of N athcrmal SAW chains of 
length L on a simple cubic lattice of M sites. The bead or 
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segment concentration is given by c = NL/M, while the 
polymer chain concentration is given by p = N/M. The 
polymers are characterized by their radius of gyration R g 
which, for an isolated polymer, scales as R g ~ L" , where 

We also define 
at which there 



v « 0.59 is the Flory exponent 
an overlap concentration p* = 1/ 



TlRl 



is on average one polymer per sphere of radius R g . So- 
lutions with p/p* < 1 are called dilute, while solutions 
with p/p* > 1 and c are called semi-dilute. When 
the monomer density c becomes appreciable, the solution 
moves from the semi-dilute to the melt regime. In this 
paper we will focus on densities p/p* < 2, i.e. the dilute 
regime and the beginning of the semi-dilute regime. 

When modelling the semi-dilute regime, it is important 
to take sufficiently long polymer chains. The first reason 
is that for studying the semi-dilute regime one needs a 
large polymer density p together with a low monomer 
density c. We found earlier that the monomer density c* 
at the overlap concentration p* scales roughly like |24|] 
c* k 4L~ - 8 for SAW polymers on a simple cubic lat- 
tice. Throughout this paper we use polymers of length 
L = 500 for which c w 0.05 at p/p* — 2, so that we 
are still clearly in the semi-dilute regime. In contrast, 
for L = 100 the monomer density is c w 0.2 at p/p* = 2, 
suggesting that a meaningful semi-dilute regime does not 
exist for such short polymers. 

The second reason for using long polymer chains is 
that we want to study properties — particularly the ef- 
fective potentials between polymer chains — in the scal- 
ing regime, where all length dependence is completely 
captured by R g . In a previous paper [ p4| we established 
that two properties of the effective potential relevant to 
thermodynamics, namely the second-virial coefficient B2 
between two polymers and the effective pair-potential be- 
tween CM's, multiplied by the square of the CM distance, 
r 2 v(r), are very close to the scaling limit for L = 500 
polymers, the length we will use in this paper. 

Simulations were done with the Monte Carlo pivot 
algorithm |^l]|3(| combined with simple translational 
moves. For concentrations p > p* we also use Configura- 
tional Bias Monte Carlo algorithms 0,|2|. For L = 500 
polymers we find that the radius of gyration of an iso- 
lated coil is R g = 16.50 ±0.03. We used a simulation box 
of size M — (240) 3 , and varied the number of polymers 
from N = 2 for the 2-body calculations to N — 6400, 
which corresponds to p/p* — 8.9. 

Since we are dealing with athermal chains consisting 
of monomers interacting only via hard-core repulsion, we 
set the reciprocal temperature [3 = 1/ksT = 1 through- 
out this paper. 



III. DENSITY INDEPENDENT MANY-BODY 
INTERACTIONS 



A. Expanding the coarse-grained free energy in a 
series of many-body interactions 

Following the discussion in |33| , the Hclmholtz free en- 
ergy T of a set of N polymers of length L with their 
centers of mass (CM) distributed according to the set of 
coordinates {r^}, in a volume V, can be written as the 
following expansion: 

N 



F{N, V, fa}) = (N, V)+J2 (r k , r ia , ) (1) 

ii < i2 

N 

+ J2 w (3) (r ll ,r l2 ,r l3 ) + ... 

il<i2<i3 
N 

... + ^ W (rn,ri 2 ...r ijv ) 



In the scaling limit, each term in the series is inde- 
pendent of L as long as the n-tuple CM coordinates 
{r^jTij ...Ti n } are expressed in units of R g , the radius 
of gyration at zero density. Note that this coarse-grained 
free energy includes an implicit statistical average over 
all the monomeric degrees of freedom for a fixed con- 
figuration fa} of the CM. The full free energy of the 
underlying polymer system can be calculated as follows: 



F(N, V) = - In J2 exp [-F(N, V, {r,})] 



(2) 



{r,} 



so that Eq. can be viewed as an expansion of the ef- 
fective interaction between the CM in terms of (entropic) 
many-body interactions. J-(°\N,V) is the so-called vol- 
ume term, the contribution to the free energy which is 
independent of the configuration {i - ;} jL4j. Here it in- 
cludes the free energy of a single isolated polymer which 
is independent of the position of its CM in a homogeneous 
solution; translational invariancc also implies that there 
is no one-body term in the expansion. Each subsequent 
term (r^ , r.; 2 . . . Tj n ) is defined as the free energy of n 
polymers with their CM positions at {r^ , Ti 2 . . . Yi n } mi- 
nus the contributions of all lower order terms. In other 
words, it is the contribution to the free energy of n poly- 
mers which is not included in the sum of all lower order 
terms. For instance, the 2-body term w^ 2 \rij) can be 
defined as the difference between the coarse-grained free 
energy T for two particles with their CM distance held 
at Tij = I r i — r j I , and the free energy of the same two 
polymers when they are infinitely far apart. Here we use 
the translational and rotational invariance of a homoge- 
neous system to reduce the number of degrees of free- 
dom. Similarly, the three-body term for a given triplet 
configuration {ri,rj,rfe} can be written in terms of only 
3 variables (see Fig. |l|) as 



,(3) 



{rij,rj kl r ki ) = T(N = 3, V, ry, r jk , r ki ) 



-J*°\N = 3, V) - w^im) - w^{r lk ) - w<V(r jk ). (3) 
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In other words, it is that part of the effective interac- 
tion between three polymers which cannot be described 
by volume and pair interaction terms alone. In princi- 
ple, this procedure may be continued until, for a system 
with N polymers, the iVth term determines the total 
coarse-grained free energy. In practice, this approach is 
not feasible because the number of n-tuple coordinates 
increases rapidly with n, as does the complexity of each 
higher order term, so that the series in Eq. (|l|) quickly 
becomes intractable. Instead, one hopes to show that 
the series converges fast enough that only a few low or- 
der terms are needed to obtain a desired accuracy. We 
now turn to the derivation of these density independent 
potentials for our system of SAW polymers. 



B. 2-body interactions 

There is a general relationship between the p — > limit 
of n-body correlation functions and the n-body potential 
|p3fl . For the 2-body case this reduces to 



lim g 2 (r) = exp 

p^O 



,(2) 



(r) 



(4) 



where g-x (r) is the pair distribution function. Although 
this definition resembles that of the potential of mean 
force (PMF), usually defined as w PMF (r) = — \n[g2(r)] 
for any density, they are only equivalent in the limit of 
infinite dilution. Strictly speaking, the PMF is not a po- 
tential but simply a restatement of the pair-correlations. 
Consider the simplest case, namely a system with no 
higher-order (n > 2) interactions. If one were to use 
a finite density PMF as a pair potential at that same fi- 
nite density, the resulting pair correlations would not be 
those of the system used to derive the PMF. In contrast, 
the potential defined in Eq. (||) is the correct pair poten- 
tial which would exactly reproduce the pair correlations 
of the original system. 

In our simulations we calculate m/ 2 -* (r) from the loga- 
rithm of the overlap probability as a function of CM dis- 
tance (c.f. Eq. |J). Although the arguments above were 
made for a free energy in a continuous space, they easily 
carry over for the lattice model we simulate. In fact, the 
CM lives on a finer grid than the original SAW lattice 
polymers, so that our results are already closer to the 
continuum limit. The overlap probability is determined 
by sampling the configurations of two polymers infinitely 
far apart with the pivot algorithm, and, after every 1000 
pivot moves, searching for any monomer overlaps as a 
function of the CM distance. The effective potential cal- 
culated in this manner has a near Gaussian shape with 
a value at full overlap of w^(0) = 1.88 ± 0.01 for our 
L = 500 polymers, very close to the scaling limit esti- 
mate w^(0) = 1.80 ± 0.05 p4L and a range of the order 

This picture agrees with 
ll studies. Note that in 



of R g p3| , ^4[|, a s shown in Fig 
earlier RG |22|] and simulation 
the scaling limit the potentials depend only on R g , so 




that the free energy cost of completely overlapping the 
CM of two polymers is independent of their length L. 
That this free energy cost at complete overlap should de- 
pend weakly on polymer length follows from their fractal 
nature [|19| , but more sophisticated scaling theory argu- 
ments ||20|| are needed to prove that w^(Q) oc L° |S,e3. 



C. 3-body interactions 

Just as the 2-body interactions follow from the low 
density limit of the pair correlations, so also the 3-body 
or triplet interactions between the CM of three polymers 
can be derived by taking the low-density limit of the 3- 
body distribution function 



33 (r 1 ,r 2 ,r 3 ) = | J ... J e -^»dr 4 .. 



.fir 



N: 



(5) 



where J-({ri}) is the coarse-grained free energy ([!]) and 
and Z = exp[— F] is the configurational integral for the 
CM coordinates defined in Eq. (||) the effective interac- 
tions already include an average over monomeric degrees 
of freedom, this configurational integral is equivalent to 
that of the full polymer system.) Taking the low-density 
limit for a homogeneous system gives 



w (3) (ri2,r 2 3,ri 3 ) = lim In 



p^O 



ga(ri2,r23,ria) 



(6) 



The (72(f) ensure that the contributions due to the pair 
interactions are subtracted from the triplet interaction 
(c.f. Eq. ||). For our homogeneous system the 3-body po- 
tential depends only on the three variables {r±2, r 2 3, ^13} 
shown in Fig. Even then, calculating the triplet in- 
teraction «/ 3 ) (ri2, r 2 3, ris) for every possible triplet ar- 
rangement is very cumbersome. We therefore confine 
ourselves to configurations that make up an equilateral 
triangle. Instead of three variables, the potential now 
depends only on the length r of each side of the triangle, 
simplifying the calculation of Eq. (|J) to: 



■/ 3) (r) = - lim 

p^O 



]ng 3 (r) - 3u> (2) (r) 



(7) 



where we also used Eq. (^). We expect that for r = 0, 
i.e. complete overlap, the three body interaction will be 
strongest, while for large r the interaction should vanish. 

We calculated u> (3) (r) for three L = 500 SAW poly- 
mers on a lattice. At this infinite dilution, 33 (r) is sim- 
ply the probability that three polymers in a configuration 
{^12, ^23, ?"3i} do not overlap. In the Monte Carlo simula- 
tion we integrate over the monomeric degrees of freedom 
by performing pivot moves. Once every 1000 MC steps 
we move the polymers into a triangular configuration and 
check for overlap. The results are plotted in Fig. ||. Since 
the total free energy increases with the number of poly- 
mers, a more relevant measure of the 3-body interactions 
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is the relative potential Au/ 3 )(r) = ui® (r)/3w^(r) 
which denotes the strength of the 3-body interaction rel- 
ative to that of the two body interactions. As shown 
in the inset of Fig. ||, the relative contribution of the 3- 
body potential to the total free energy is quite small, only 
about 9% of the contribution from the pair potentials. 



D. 4-body interactions 

Next we turn to the 4-body interactions. Again, even 
for a homogeneous phase, the total number of relevant 
coordinates makes the calculation of the full interaction 
prohibitively complex. To restrict the number of coor- 
dinates in our calculations we determine the 4-body po- 
tential by placing the four polymer CM's on a regular 
tetrahedron and determining the non-overlap probability 
as a function of the length r of each side of the tetrahe- 
dron. The 4-body potential is then defined as 

w (i) (r) = - lim In \g 4 (r)} - 4w (3) (r) - 6w (2) (r) (8) 

since the tetrahedral configuration includes four equilat- 
eral triangles (3-body interactions) and six edges (2-body 
interactions). The full and relative 4-body interaction, 
Aw (4) (r) = u> (4) (r)/(4w (3) (Y) + 6u; (2) (r)), are plotted in 
Fig. ^. Note that the 4-body interaction is smaller in ab- 
solute magnitude and has the opposite sign to the 3-body 
interaction. The relative contribution of the 4-body in- 
teractions to the total free energy is less than 5% of the 
total potential, and also less than the relative 3-body 
contribution. 



E. 5-body interactions 

Calculating the 5-body interaction is even more com- 
plicated than the 4-body interaction, and so we only eval- 
uate it at full overlap - when all the CM's coincide - 
where we expect its contribution to be largest. More 
generally, for any nth order term the interaction at full 
overlap is given by: 



,(n) 



(0) 



lim[ln ff „(0)] 



m—n — l 

E 

m=2 



w (m \0) (9) 



where lim^o SVi (0) is the normalized probability of full 
overlap of the CM of n polymers. As long as the par- 
ticles are equidistant from each other, the same combi- 
natorial expression holds for finite r. For the 5-body 
term we find that u/ 5 )(0) = —0.4 ± 0.15, while the 
relative contribution of the 5-body terms is given by 
w^{0)/{5w^{0) + 10w (3 \0) + 10w^(0)) = 0.027±0.01. 
Again, the relative contribution of the 5-body term to 
the free energy is smaller and of opposite sign to those 
of the 4-body terms. Going beyond the 5-body inter- 
action, even at complete overlap becomes increasingly 



difficult. For example, for the 5-body interaction, of 
the 10 s overlap checks, each attempted after 1000 pivot 
moves, only about 30 resulted in non-overlap. For the 
6-body interaction we estimate that 10 11 MC attempts 
would be needed. Another problem arises from finite 
monomer density. As more and more polymers overlap, 
the monomer density increases, so that in practice for a 
given polymer length L, only a finite number of multi- 
ple overlaps are possible. We found previously that the 
largest finite-size corrections to the scaling limit were at 
r = for u)W(0) pi. The same probably holds for the 
higher order interactions. However, our L = 500 calcula- 
tions should still be very near the scaling limit. 



IV. EFFECTIVE DENSITY DEPENDENT PAIR 
INTERACTIONS 

From the previous section we see that explic- 
itly calculating the density independent interactions 
w 



in) 



(r ii; r.j 2 . . . v iN ) becomes rapidly more complex with 
increasing order n. Calculating all higher order terms is 
therefore impossible. In this section, we describe a way 
to include the average effect of all higher order terms by 
extending the relationship between the pair-interactions 
and the pair-correlations to finite density p. This leads 
to a density dependent effective pair interaction v(r; p). 



A. Inverting pair-correlations to derive density 
dependent pair potentials 

Although for finite densities there is no known direct 
functional relationship of the type of Eq. (^), there is a 
theorem which states that for any given pair-correlation 
function g<i (r) and density p, there exists (up to an ad- 
ditive constant) a unique pair potential v(r; p) which ex- 
actly reproduces 172(f) at that density Q,^5|. If the orig- 
inal (72(f) is generated by a system with only pair inter- 
actions, then v(r; p) = ur 2 )(r) will be independent of 
density. If there are any higher order interactions in the 
original system which influence the structure, then this 
equivalence will only hold for lim^o v(r; p) = w^ 2 \r). 
At finite densities v(r; p) must change since the struc- 
ture is no longer equal to the one generated by m/ 2 -* (r) 
alone. Therefore v(r; p) must be density dependent. 

In fact, this is what we found in two previous pa- 
pers |3^,|2J], where we used the hypernetted-chain (HNC) 
approximation from liquid state theory [3^ ] to extract 
v{r; p) from computer simulations of the <?2 M's between 
the CM's of SAW polymers. For completeness, we show 
these effective pair interactions in Fig. (||). As expected, 
there is a clear density dependence. Without going into 
much detail about the inversion of v(r; p) from (72(f), we 
do want to point out that the process can be very subtle. 
As illustrated in Fig. |], the 52 Ws generated at p — p* by 
v(r; p — 0) and v(r; p = p*) are very similar. Any tech- 
nique to derive v(r; p) from 172(f) must be significantly 
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more accurate than the difference between the g2(r)'s 
shown in the figure. The accuracy of the techniques we 
use has been discussed in Ref. [M , and will be analyzed 
in much more detail in another publication |]37f . 

Any approximation which correctly reproduces the 
pair-correlations will also predict the correct thermody- 
namics through the compressibility equation |36|,[58| . For 
our density dependent v(r; p) this is indeed the case, since 
we found good agreement between the equation of state 
(EOS) Il/p generated by the effective potentials in Fig. || 
and the EOS of the underlying SAW polymer solution. 
In contrast, the v(r\ p = 0) potential underestimates 
the EOS, and we find mean- field fluid behavior [||J40| 
II/p ~ p at large p instead of the correct 11/ p ~ /r 73 
scaling. So, even though the p = potential results in 
pair-correlations g 2 {f) that are similar to the true g 2 (r)'s, 
the effective thermodynamics can differ significantly. The 
difference arises from the neglected many-body interac- 
tions, as discussed in the Appendix. 



B. Understanding the density dependence of the 
effective pair-potential 

Given the success of the density-dependent pair inter- 
action in describing pair-correlations and thermodynam- 
ics, we next turn to the question of whether the density- 
dependence of v(r;p) can be directly understood from 
the density independent many-body interactions. 

Within the HNC approximation, the following expres- 
sion due to Reatto and Tau || and also Attard || 

v(r 12 ;p) = w^(r 12 )-p J (e-^ 3) (n2,r 13 ,r 23 ) _ ^ (1Q) 

xg2(r 13 ; p)g 2 (r 23 ; p)dr 3 , 

describes the density dependence of the pair-potential 
that would reproduce the true pair-correlations induced 
by the 2-body and 3-body potentials. This is a general- 
ization of earlier expressions |4l]j4^] and neglects terms of 
order 0(p 2 ) and higher. In the literature it has mainly 
been applied to the Axilrod- Teller interaction for rare- 
gas fluids, where it works remarkably well, see e.g. |43j | 
and references therein. 

Figure || highlights the density dependence by plot- 
ting i> ex (r; p) = (v(r; p) — w^(r))/p. For clarity, we 
have replaced the rather noisy data by spline fits. For 
densities p/p* < 1 the curves are close to each other 
suggesting that the roughly linear density dependence in 
Eq. ([l0|) holds true. For larger densities into the semi- 
dilute regime, v ex (r] p) becomes smaller in magnitude 
and the density dependence becomes non-linear. This 
non-linearity is not unexpected, since Eq. ( |l0| ) neglects 
higher order terms in p, as well as the effects of 4-body 
and higher order interaction terms. 

We can go even further and directly calculate the 
triplet induced density dependent term v ex (r;p) by sub- 
stituting Eq (||) into Eq. ( |i~0| ) to obtain 



v ex (ri 2 ,p) 



lim 



33(^12, n 3 ,r 23 ) 



v p^o g 2 (ri 2 )g 2 (n 3 )g 2 (r 23 ] 
xg 2 (r 13 ; p)g 2 (r 23 ; p)dr 3 . 



- 1 



(11) 



The g 2 (r]p), in contrast to the g 2 (r) in the first term 
in (pd|), are defined at the density of interest. Evalu- 
ating this integral is difficult, because the term between 
brackets can become very small. We use a direct MC pro- 
cedure, where two polymer coils are held with their CM 
a distance r\ 2 apart while we integrate over the position 
of the third particle. In order to ensure that the inte- 
gral converges it is crucial to use the g 2 (r) at p = (i.e. 
those between the brackets in Eq. |ll|), from the simula- 
tion itself, by calculating it on the fly. This is necessary 
to avoid small errors in the radial distribution function 
which build up during the integration over the volume. 
The g 2 (r; p) at finite density are known from previous 
calculations. The resulting v ex (r; p) for p = is plotted 
in Fig. ||. The results for finite density do not differ by 
very much. 

In conclusion, the density dependence is mainly caused 
by 3-body interactions, at least in the dilute regime. At 
semi-dilute densities, higher order many-body interac- 
tions may come into play. 



V. MANY-BODY CORRELATIONS 

If g 2 (r) is generated by a system with only pair- 
potentials, then the exact inversion of g 2 (r) at any den- 
sity will reproduce the exact pair potential. For such 
a system, the inverted pair potential can be used in 
principle to determine all higher order correlation func- 
tions of the original system. In other words, for systems 
with only pair-interactions, the pair-distribution function 
g 2 (r) contains enough information to generate all higher 
order correlation functions [Ql 

If a system has 3-body or higher order interac- 
tions, then our v(r; p) still exactly reproduces the pair- 
correlations. But, as we shall demonstrate in this section, 
it can no longer exactly reproduce the higher order cor- 
relations. Nevertheless, we will show that the differences 
are not very large for the case of SAW polymers at the 
densities we study. 



A. Bond angle distribution from 3-body correlations 



Calculating and comparing the full 3-body correlation 
functions would be very cumbersome for many of the 
same reasons that it is difficult to map out the full 3-body 
interaction. Therefore, we resort to a reduced picture 
where a subset of the variables are integrated out |45| , fl6f . 
One popular measure of the 3-body interactions is the 
bond angle distribution function, defined as 



G 



b(0, r c ) = 8ir 2 p 2 N f ' f g 3 (r 12 ,r 13l (rf 2 + r 2 3 (12) 
Jo Jo 

— 2ri2ri 3 cos9) 1 / 2 )rl 2 rl 3 sin9dri2dri3, 

where N is a dimensionless normalization constant. This 
integral sums over all triplets within a cutoff radius r c 
from the central particle and determines the distribution 
of the angle 9 in these triplets. We calculated the bond 
angle distribution for both the SAW simulations and the 
effective pair potentials for different cutoff radii r c as 
shown in Fig. ^. The effective potentials v(r; p) reproduce 
this measure of the 3-body correlations remarkably well. 
Since for an ideal gas the bond angle distribution exactly 
follows a sine curve, dividing the bond angle distribution 
by sin# highlights the deviations from ideal behavior. In 
Fig. ^ we show the renormalized bond angle distribu- 
tion. The differences between the curves are now clearer. 
The absolute deviations from the sine-like behavior are 
largest at small 9 because the particles repel each other 
and triplets with small 9 will be relatively rare. In the 
case of a hard sphere systems this correlation hole would 
be even more pronounced. It is remarkable how well the 
bond order distribution follows the ideal sine curve for 
the larger angles. At a cutoff radius of r c = 0.5 the dis- 
tribution becomes flatter, reflecting the broad flat top of 
the repulsive Gaussian shaped pair potential. 

Instead of determining the distributions from explicit 
simulations, we can also substitute the Kirkwood super- 
position approximation pq 



ff3(«"i.2,ri3,r23) ~ 92^12)92^13)92^23) 



(13) 



into Eq. [12] and calculate the integral directly by using 
the radial distribution functions from previous simula- 
tions. This approximation is also included in Fig. ^| and 
Fig. and turns out to be very accurate, except for 9 m 
and 9 w ir (see Fig. ^) where the simulations are prone 
to large statistical errors, due to the vanishing volume of 
the available phase-space. 

The bond-angle distribution is not very sensitive to dif- 
ferences in the full 3-body correlations (see e.g. Ref. [Kjfl), 
partially because it is an integrated quantity. An ex- 
ample of this is given in Fig. || where the bond angle 
distributions of the effective potentials v(r; p = 0) and 
v(r;p = p*) are compared for the same density p — p* . 
Clearly, there is hardly any difference between the distri- 
butions. 



B. Deviations from Kirkwood superposition for 
3-body correlations 

A more sensitive measure of triplet correlations is the 
deviation from the Kirkwood superposition approxima- 
tion, Eq. (113), which we define as: 



Since our v(r; p) exactly reproduces the (feW, this ex- 
pression should highlight any differences between the true 
g 3 and the g 3 arising from our effective potential picture. 
To simplify, we limit ourselves for a given r%2 to triplet 
configurations for which r = r'23 = r± 3 (i.e. isocclcs tri- 
angles) . 

First, we compare in Fig || the G 3 (r) at p = p* gen- 
erated by v(r;p — 0) and by v(r; p — p*). Just as we 
found for the bond angle distributions, the G 3 (r) are very 
similar even though the potentials are different. We al- 
ready showed that the 52 (?") are not very different either 
(see Fig. ||) , so that the same now holds for the full 
g3(ri2,r 23 ,r 13 ). 

Next, we turn to a comparison between the true 
Gf AW (r) derived from explicit simulations of our SAW 
polymer system and the G§ ff (r) of the effective potentials 
at p = p* . As can be seen in Fig. ||, our effective pair- 
potential v(r; p) does not exactly reproduce the SAW 3- 
body correlation function. The trends are similar, but 
the deviation from superposition of the SAW polymers, 
Gf AW (r), is consistently larger than the same quantity 
generated by the effective potentials G 3 s (r), especially if 
r < R g . 

For systems with an explicit three-body interaction 
the Kirkwood superposition approximation is sometimes 
written as: 

93(r '12, r '13, r 23 ) « 02 (»"12)52 (13)02 (ras) (15) 
xexp -w (3) (ri2,ri3,r 23 ) . 

This is exact in the p — » limit, as can be seen 
from Eq. (||). Note that in this same limit the 
three-body correlations induced by the v(r; p) reduce 
to the simpler Kirkwood superposition approximation 
53(^12, ^13, ^23) ~ g2{ri2)g2{ri 3 )g2{r2 3 ) 1 demonstrating 
explicitly that in contrast to the 2-body correlations, 
v(r; p) cannot exactly reproduce the 3-body correlations 
if there is a 3-body interaction present! In fact, we have 
shown explicitly that two systems with identical pair cor- 
relations, namely our effective potential system and the 
original SAW system, can have differing triplet correla- 
tions. 

These arguments also suggest that a simple approx- 
imation, namely G 3 AW (r) w G 3 S (r) exp[— (r)], can 
shed some light on the differences observed in Fig ||. 
Since (r) is negative for equilateral triangle config- 
urations, as illustrated in Fig|^, it is perhaps not surpris- 

for the 



eff 



G 3 (ri2,r 13 ,r 23 ) 



gz{ri2,ri 3 ,r 23 ) 
92{ri2)g2(ri 3 )g2(r2 3 ) ' 



(14) 



ing that roughly speaking G§ AW (r) > G 
isoceles triangle configurations plotted in Figs |9|, at least 
in the region where (r) is non-zero. Unfortunately 
the statistical errors in this region are very large, making 
a quantitative comparison difficult, but the deviation is 
certainly of the same order as would be expected from 
an extra factor exp[— (r)] (compare with Fig.||). 

Our very simple approximation illustrates how devi- 
ations from Kirkwood superposition originate from two 
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effects: (1) Deviations induced by correlations generated 
by the pair-potentials alone. These have been studied in 
great detail for hard-sphere systems, see e.g. (2) 
Deviations induced primarily by 3-body potentials. In 
practice, of course, these two effects are somewhat en- 
twined, especially at higher densities. Nevertheless split- 
ting the two effects can shed light on the origin of three- 
body correlations. In particular, it suggests that while 
an effective pair-potential v(r; p) that exactly reproduces 
the #2 (r) can partially reproduce deviations from super- 
position of type (1), it will fail for deviations of type 
(2). 

Since the 3-body and higher order interactions between 
the CM of polymer solutions are not very strong, the to- 
tal <?3(r)'s are still remarkably well reproduced by the 
v(r; p), especially when integrated quantities such as the 
bond-angle correlations are considered. However, for sys- 
tems where 3-body interactions are strong, such as liquid 
Si or liquid Ga, one cannot expect the same success from 
effective pair potentials. Very similar conclusions were 
stressed by Evans B] in the context of reverse Monte- 
Carlo simulations EcTIiJl. 



VI. CONCLUSIONS 

Integrating out the monomeric degrees of freedom to 
obtain a description based on effective potentials between 
polymer CM's is a useful coarse-graining technique for 
polymer solutions. Because simulations can be performed 
to high accuracy, the lessons learned here should be appli- 
cable to a much broader range of coarse-graining schemes. 

In particular, we showed that the free energy of the 
polymers can be expanded in a series of state indepen- 
dent many-body effective potentials. The terms in the 
series oscillate as (—1)", and become smaller in absolute 
magnitude for increasing n. The scaling theory devel- 
oped in the Appendix confirms these ideas, and can be 
used to extend them to arbitrary order n. 

A parallel description of the coarse-grained polymer 
solution was developed in terms of an effective state (den- 
sity) dependent pair potential v(r;p), which exactly re- 
produces the pair-correlations and, in an average way, in- 
cludes all the higher order terms in the many-body free 
energy expansion. The density-dependence of this effec- 
tive pair-potential can be largely understood from the 
direct influence of the density-independent three-body 
interactions. 

The three-body correlations are also well described 
by this effective pair-potential picture. If the bond an- 
gle distribution is used, the differences between the full 
SAW polymer system and our effective pair-potential pic- 
ture are almost indistinguishable. However, this is not 
a good "order parameter" for measuring deviations in 
three-body distribution functions, since there is almost 
no difference between results from the full simulations 
and those produced with the Kirkwood superposition ap- 



proximation. When we use a more direct measure of 
the deviations from Kirkwood superposition, small differ- 
ences between the effective two-body and the full SAW 
triplet distributions can be measured. These arise pri- 
marily from the direct effect of the three-body potential 



IV 



(3) 



(r), and illustrate a more general point, namely that 
an effective two-body interaction can never simultane- 
ously reproduce the two and three-body correlations in 
a system with many-body interactions exactly. Since for 
polymers in a good solvent these many-body interactions 
are relatively weak, a coarse-grained description based on 
effective pair-interactions works remarkably well, at least 
for the dilute and the beginning of the semi-dilute regime 
p3| , pi[ . Whether this success can be extended deeper into 
the semi-dilute or into the melt-regime, or even for poly- 
mers in poor solvents, remains to be seen, and will be the 
subject of future investigations. 
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APPENDIX A: RELATIVE STRENGTH OF N 
BODY TERMS FROM SCALING THEORY 



In a recent paper, von Ferber et at 25 used scaling 
theory and simulations to calculate the triplet interac- 
tion for star-polymers and found an attractive interac- 
tion with a relative strength of about 11%, very similar 
to our 3-body results. The natural choice for the po- 
sition coordinate of a star-polymer is not the CM, but 
its midpoint |jo|. However, we will see that for esti- 
mating relative contributions, the difference between our 
CM representation and the midpoint representation is 
not that important. 

Here, we apply the star-polymer scaling theory to esti- 
mate the relative contributions of w n (0) to all orders in 
n. We specialize to linear polymers, which can be seen as 
star-polymers with only two arms. We first note that the 
partition function for n polymers with their mid-points 
constrained to be at a distance r <c R„ apart scales as 



5p 52 1 



Z n (r) 



(Al) 



in the limit r/R g — » 0. Here 9 n is the contact exponent 
which in turn can be written as 



= nrj 2 - 772- 



(A2) 
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where the rjf are the scaling exponents for a star-polymer 
with / arms. These are tabulated for two different renor- 
malization group calculations in [^5],[33| , and can also be 
approximated by a simpler expression 



n f w 0.335307/ 



.3/2 



(A3) 



which is expected to become more accurate for larger / 
|54j| . Comparing the different approximations gives an 
indication of the accuracy of the scaling theory results. 

The probability of finding n polymers with their mid- 
points a distance r <C R g apart can be found from the 
partition functions since: 



lim ln< id (r) oc In [Z n {r)\ 

p^O 



Mr/Rg 



(A4) 



Here g™ td (r) is the n-body distribution function for the 
midpoints of n polymers. These scaling theories can be 
developed for polymers in the midpoint representation in 
part because of the analogy of the overlap of n polymers 
to a In arm star-polymer. Such an analogy is not as ob- 
vious for the CM representation, hampering the deriva- 
tion of a CM based scaling theory. On the other hand, 
MC simulations for overlap in the CM representation are 
much easier than for the midpoint representation because 
the probability g n (r) is much larger than g™ %d {r) in the 
small r limit, in fact g™ ld (0) = while the CM g n (0) is 
finite. Gathering good statistics for the mid-point repre- 
sentation is thus much slower than for the CM represen- 
tation. It nevertheless seems reasonable to expect that 
the relative strengths of different orders of the interac- 
tions in the CM representation are similar to those of 
the midpoint representation, suggesting that the relative 
probability of full overlap of the CM of n polymers scales 
as: 



m<?n(0) 



lim . 

p^o \\ng m {0) 



lim lim 

r-K) p^O 



in gr d (o 

\ng™ id (r) 



2- (A5) 



In Table [j] we confirm this ansatz by comparing 
lng„(0)/lnf/2(0) calculated by direct simulations for n — 
3, 4, 5 with # n /#2 calculated by three different versions of 
the scaling theory. Clearly, all three scaling theories and 
the simulations agree reasonably well with each other, 
giving us confidence to proceed. 

Armed with this approximate equivalence of the g n 's 
in the two representations, we can now fruitfully com- 
pare the strength of the n body interaction w^ n \0) for 
CM particles to an expression derived from the scaling 
theory, by using Eq. (A2) and Eq. (A4) to recursively 
rewrite Eq. (0) as: 



wW(o) oci&M = (-i)« {^j(- l ) m ^m-m- m ). 



m=2 



(A6) 



Here is the coefficient of the full (divergent) 

midpoint-midpoint interaction —w^h(r) ~ w^\a{r/Rg) 



in the limit r — + 0. The values of these terms relative to 
the n = 2 term are compared in Table || with direct simu- 
lation results for n — 3, 4, 5, and again good agreement is 
obtained. For the scaling theory each higher order term 
is opposite in sign and smaller than the previous one. 
Our simulations show the same trend for the sign, but 
the error bars on the n — 5 simulations are still too large 
to ascertain that its magnitude is less than the n = 4 
term. 

The strength of the n body term contribution to the 
free energy relative to all lower order terms can also be 
found from the arguments above. In this case the pre- 
f actors cancel and we find: 



Au/™)(0) « At&W 



(A7) 



Once again the simulations and the scaling theory agree 
quite well both for the magnitude and the sign of the 
different terms. 



APPENDIX B: CONVERGENCE OF THE FULL 
FREE ENERGY SERIES 

The good agreement found between the scaling the- 
ory and the simulations for n = 3, 4, 5 suggests that we 
can evaluate trends for arbitrary order n from the scaling 
theory. We use the simple expression in Eq. (A3) which 
gives r/f for arbitrary /, to simplify Eq. (A6) to 



iiW = 0.9484(-l) n (J>2 ^j(-l) m+1 m 3 ^ • ( B1 ) 

The sum between brackets is positive for all n, and be- 
comes smaller for each subsequent term. For example for 
n = 3 it is 0.29, while for n = 10, 000, it is has dropped 
to 0.009. This implies that at all orders, the n-body 
interaction w^ n '(Q) oscillates in sign as (—1)" and de- 
creases in absolute magnitude with increasing n. The ab- 
solute magnitudes of the relative contributions Abj'"' (0) 
decrease even faster with n. 

We expect that the absolute magnitude of each n body 
term is largest at full overlap. However, the fact that 
each higher order term decreases in absolute magnitude, 
or even that the terms oscillate in sign at r = 0, does 
not necessarily imply that the full free energy expan- 
sion of Eq. (|l|) will converge quickly for all configura- 
tions {rj} of the CM. For one thing, for an N particle 
system, the number of n-tuple coordinates grows as the 
binomial N\/(N—n)\n\. Therefore, even though the mag- 
nitude of the n body terms (r^ , rj 2 . . . ) becomes 
smaller with increasing n, the number of such terms at 
each order n in the free energy series Eq. ([l]) increases 
almost exponentially with increasing n. Of course for 
a given configuration {r^}, not only the magnitude, but 
also the range of each term w^ n ' (r^ , r^ 2 . . . rj n ) is needed 
to decide how many, if not all iV terms, are needed to 







calculate the full free energy to a desired accuracy. As 
mentioned before, calculating the complete dependence 
of w^ n ' (rjj , r,- 2 ...Ti n ) on all possible n-tuple CM coor- 
dinates {r^ , Ti 2 . . . Ti n } is usually an impossible task for 
higher order n. 

Nevertheless, some progress can be made by looking 
at the special case of a set {r^} where all N CM's are at 
the same point. This is simply the situation studied in 
the previous section by scaling theory, for which the free 
energy is given by 



T(N, V, {n}) - J*°\N, V) = - lng N (0) oc 



'JV< 



(B2) 



Ignoring for the moment the volume term, who's con- 
tribution is expected to be negligible Jlj],|24| we can ask 
the question: How well is the free energy described by 
taking only pair-interactions terms into account? Since 
there are ^N(N — 1) pair terms, the free energy with 
only pair-terms taken into account scales as T oc N(N — 
l)u>( 2 '(0) ~ iV 2 while the true free energy scales as 
J~ oc 9n ~ TV 3 / 2 . Including only the pair terms heav- 
ily overestimates the free energy, a result also found by 
von Ferber et al. |25f| . In other words, if one truncates the 
coarse-grained free energy series in Eq. (Q) then certain 
configurations {r^} will not be properly accounted for. 

However, if one were to use this coarse-grained free en- 
ergy series to study a given polymer solution, then the 
relative probability for finding a configuration at com- 
plete overlap of the CM would be very small. On the 
other hand, for any given order n, the probability of find- 
ing n-fold overlap increases with increasing density, so 
that one would expect to need more and more terms in 
the coarse-grained free energy series as the density in- 
creases. 

Another way to measure the effects of truncating the 
coarse-grained free energy series in Eq. ([lj) would be to 
compare it to the excess free energy per particle for the 
original polymer solution. This can be found by averag- 
ing T{N, V, {i"i}) over all configurations {riito obtain 
the total free energy F(N, V), as done in Eq.(J2J). Includ- 
ing only pair interactions in our CM representation of 
polymer solutions casts the problem into that of finding 
the free energy of a mean-field fluid [ p9| , p3|j4C| ] for which 
the excess free energy per polymer scales as: 



pex 



(B3) 



while the true excess free energy of a polymer solution in 
the semi-dilute regime follows from scaling theory Eg] 



pex 

~N~ 



,1-3 



(B4) 



Whereas including only the pair interactions in an aver- 
age over all configurations {r^} in the semi-dilute regime 
underestimates the total free energy, taking into account 
only pair terms for single configuration {r,} at com- 
plete overlap of the CM overestimates the coarse-grained 



free energy. This implies that the special configurations 
where the coarse-grained free energy series in Eq. (|l|) 
breaks down do not have a strong influence on the ther- 
modynamics. 
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TABLE I. Comparison of scaling theory and simulations for many-body interactions between the polymer CM. The labels 
a and b denote results that follow from two different renormalization group calculations for the exponents rj n Pq| , while label 
c comes from the simple expansion of Eq.b^3J, and label d denotes simulation results for L — 500 SAW polymer simulations. 



R 


2 


3 


4 


5 


6 


7 8 9 10 


a 9 n /8 2 


1(9 2 = 0.8) 


2.55 


4.49 


6.75 


9.3 




b 


1 (ft, = 82) 


2.65 


4.80 


7.41 


10.43 






1 (0 9 = 791 


2.65 


4.83 


7.46 


10.50 


13.90 17.66 21.73 26.10 


d In Or, (0) 1 In ao(0) 


1 (lno?(0) = 1 88 


) 2.74 ± 0.01 


5.13 ± 0.02 


7.99 ± 0.05 






a w (n) /w (2) 


1 (w {2) = 0.8) 


-0.45 


0.29 


-0.19 


0.11 




b 


1 (i& (2) = 0.82) 


-0.35 


0.22 


-0.15 


0.085 




c 


1 (w (2) = 0.79) 


-0.35 


0.22 


-0.17 


0.14 


-0.12 0.11 -0.10 0.094 


d w (n) (0)/w (2) (0) 


1 (w (2) (0) = 1.88 


) -0.26 ± 0.013 


0.17 ± 0.031 


-0.21 ± 0.08 






a Am"(0) 




-0.15 


0.069 


-0.027 


0.012 




b 




-0.12 


0.048 


-0.019 


0.0083 




c 




-0.12?? 


0.049 


-0.023 


0.014 


-0.0088 0.0063 -0.0046 0.0036 


fl 




-0.09 ± 0.004 


0.034 ± 0.007 


-0.026 ± 0.01 
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FIG. 1. The three variables {7*12, ri3, 7-23} which char- 
acterize a triplet configuration. The angle 9 between 
the vectors ri2 and ri3 is related to the distance r23 by 

»"23 = r 12 + Tl 3 " 2ri2ri3 COS0. 





1 2 3 4 5 

r/R 9 

FIG. 3. The effective polymer pair potentials v(r;p) de- 
rived at different densities from an HNC inversion of the CM 
pair distribution functions g%(r) of polymer coils depicted in 
the inset, (from Ref @) 




FIG. 2. Effective potentials 7/j (n) (r) for n = 2 (solid line), 
n = 3 (dashed line), n = 4 (dash-dotted line), and n = 5 
(symbol with error bar at r = 0). Inset: relative potentials 
Aw' n '(r) for n = 3 (dashed line), n = 4 (dash-dotted line), 
and n = 5 (symbol with error bar at r = 0). 
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=0.81 


p/p- 


=1.09 


p/p> 


=1.64 


p/p> 


=2.18 



-1.5 



3 

r/R. 



FIG. 5. The density dependent excess potentials v ex (r; p) 
determined for several densities from the SAW simulations 
are denoted by the smooth spline fit curves. The open circles 
with error bars correspond to the evaluation of Eq. ( |ll| ) at 
p = 0, the line through them is to guide the eye. 




0.2 0.4 0.6 o 

6/71 



0.2 0.4 0.6 0.8 1 
8/7C 



FIG. 7. The normalized bond angle distribution b{9, r c ) 
plotted for several cutoff radii r c for p = p* . The solid curves 
denote the effective potential results, the dotted curves corre- 
spond to the explicit SAW simulations and the dashed curves 
show Kirkwood's superposition approximation. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

2 | . 1 . 1 . 1 . 1 . 1 2 i . 1 . 1 . 1 . 1 . 1 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



6/71 8/71 

FIG. 6. The bond angle distribution b(6, r c ) plotted for 
several cutoff radii r c for p = p* . The solid curves denote 
the effective potential results, the dotted curves correspond 
to the explicit SAW simulations and the dashed curves show 
Kirkwood's superposition approximation. 




0.2 0.4 0.6 0.8 1 



e/7t 

FIG. 8. Comparison between bond angle distribution 
b(9,r c ) for effective potentials v(r\ p == 0) (dotted lines) and 
v(r; p = p*) (solid lines), both determined from simulations at 
p — p* . From left to right the curves correspond to r c — 2.5, 
r c = 2.0 and r c — 1.5. The differences are so small that the 
lines are virtually indistinguishable on the scale of the graph. 
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FIG. 9. To measure G2,(r\2, ?"i3, r\s), which denotes deviations from Kirkwood superposition, we fix ri2, and plot the Gs(r) 
as functions of r = ri3 = r23 (isoceles triangles, see Fig. |l]). Results are given for SAW polymers (circles), effective potentials 
v(r; p) (solid lines), and v(r; p — 0) potentials (dashed lines). All plots are from simulations at p = p* . Note that the two 
effective potential plots are more or less identical within the statistical noise, while the SAW G§ AW (r) is slightly higher by 



approximately a factor exp(- 



,0) 
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